1 



Effective Macroscopic Response of a Composite with Small Deviations 
from Periodicity: Application to Colloidal Crystals 

Sergey V. Barabash, David Stroud a 

a Department of Physics, The Ohio State University, Columbus, Ohio 43210 



Using the spectral approach, we analyze the effective properties of a composite which deviates slightly from 
periodicity. We find that, when the inclusions are randomly displaced from their equilibrium positions, the sharp 
resonances seen in the periodic case are broadened, and an additional branch cut appears. We use these results 
to analyze the effective dielectric constant of a colloidal crystal. 



The spectral approach[[l| has proven useful in 
analyzing the effective dielectric constant of a 
composite material. For example, it has led to 
the derivation of exact bounds on such effective 
properties [^H, as well as to approximate calcu- 
lations of effective nonlinear properties H; and it 
can be generalized to polycrystalline materials 
However, it has proven difficult to apply to the 
spectra of real systems, with the exception of a 
few simple cases, such as a simple cubic lattice of 
identical spheres 

In this paper, we achieve some progress in this 
direction by calculating the effect of weak disor- 
der on a periodic arrangement of identical inclu- 
sions. Our results may prove useful in describing 
the effective properties of real systems such as 
colloidal crystals. 

First we review the periodic case. Consider a 
system of identical inclusions (labeled by the in- 
dex a). Let the spectrum of each inclusion be 
described by the poles {s aa } — {s a } with the 
corresponding amplitudes {M aa } — {Al a }, and 
let Q aa ,/3b be the overlap integrals between the 
inclusions (see Ref. jjj for further definitions). 
The poles of the system are then given by the 



solutions of 

(s — S aa )A aa = ^ Qaa,f3bA{3b, 
0,b^a 

where the amplitude of the pole s^is 



(1) 



(2) 



In the purely periodic case, Bloch theorem implies 
that A aa = AQ,(k)e lk ' Ra , where R a is the position 
of the inclusion a. Since Q aa , pb = Qa^R-b — R-a), 
it then follows |l|] that the spectrum is a system 
of "bands" given by 



(s - s a )A a {k) = Qa,/3(k)^(k), 



where 

Q af3 (k) = Y^Q^pbe'^ ■( R *- R »>. 

b^a 



(3) 



(4) 



For the periodic case, the entire amplitude within 
each band is concentrated in the k = pole: 



M«(k) 



(5) 
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where N — > oo is the number of inclusions in the 
system. 

We now introduce disorder into the system by 
displacing the inclusions to positions 



R„ 



ikR a 



(6) 



Eq. (^) could, for example, describe a system of 
phonons in a colloidal crystal. Each Uk defines 
the amplitude and polarization of a phonon with 
wave vector k. We will assume that all nonzero 
Uk's are of the same order of magnitude: |Uk| ~ 
U p h, and treat U p h as a small parameter. To find 
the effect of such a "phononic" perturbation on 
the spectrum of the system, we express A aa as 



A a a — ^ A a (k)e ** , 



(7) 



where k is confined to the first Brillouin zone. 
Substituting (0) into ([!]), we find that Eqs. | 
should now be generalized to 



(s - s a )A a (k) = J2 Q*A k i k'M/3(k') 



(8) 



/3k' 



Q^(k,k') = I^^Q aai/36e Kk' R b -k.R a ) (9) 

a b=jta 

If we restrict our consideration to the correc- 
tion of the lowest order in U p h, we can explicitly 
write down the approximate expression for the 
position-dependent overlap integrals Q aa ,pb'- 

Qaa.fib = Qa,/3(R-b — R-a) + 

80 {1) 



fib' 



U k ' (VQ a , /3 (R b -Ra)) 



k R h 



ik-R„ 



) (1 \b 
" aa,po 



(10) 
back into 



Plugging in this expression for 
Eq. (H), we obtain: 

g Q ^(k,k') = S Kk ,Q a p(k)+ U k -k"Qa,/3(k,k'),(ll) 

where 

Q Qi/3 (k, k') = (VQ^) k - <yQ a ,p)» (12) 
and 

(VQ a ,/j)k= E [VQ Q ,MAR)] x e ikAR , (13) 

AR^O 



with AR = Rb — R a - Note that because 
the Fourier transforms ([|) and ( |l3|) involve only 
discrete sums rather than continuous integrals, 
one cannot assume that (V(5 a A = kQ a ,p(k). 
Instead, both (Q) and (|l3| ) should be evalu- 
ated directly from the known functional form of 
Q a ,/j(AR). 

The Eq.(||) now takes the form 



(s - s a )A a (k) = |Qa, / 3(k)A (3 (k) 
+ ^ u k-k'-Qa,/3(k,k')A /3 (k')} 



(14) 



The solutions of this equation give the new po- 
sitions of the poles and the corresponding eigen- 
vectors. Note that unlike Eq. (^), which could be 
solved for each value of k independently, a sin- 
gle solution of ( |l4| ) generally involves an infinite 
number of amplitudes A a (k'). 

We now explicitly find the solutions of ( |l4| ) , re- 
stricting ourselves to the case when in the purely 
periodic system the pole at Sq (ko) is not degen- 
erate. We can expand the corresponding eigen- 
vector Aa' k °\k) of Eq.(14) in terms of the un- 
perturbed eigenvectors Aq-J^ (k) by writing 

4p»)(k) = Ao«(k)^k +E4!k ko) ^oi j) (k),(15) 



where the 



cj I k k °' ) are of the first order 
in phononic amplitudes. This can be done be- 
cause(|] Q* a p{k) — Qp^ a (k), and thus the dif- 
ferent solutions Ao^\k) of (||) do form a com- 
plete set at each given value of k; we also assume 
that J2 a ^oL 4) (k)*^oL j) (k) = S id Using the latter 
property and noting that 







)/3 (k) - s a 5 a ,/3 A Q ( ^(k) 



= 4 4) (k)A «(k) (16) 

[cf. eq. (^)] it is straightforward to show that the 
first order correction to the eigenvector is given 
by 



„(i,ko) 
C i,k 



(17) 



a, 



Q Q ,Mk,k )4^(k ) 



S «(ko)-4 J) (k) 
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for k 7^ k . For k = k , Cj^' = 0, reflecting the 
fact that the diagonal term of the perturbation 
is zero: Q Q , j( g(k, k) = 0. For the same reason 

the position of the non-degenerate pole s (ko) 
remains unchanged in the first order. 

Next, we analyze how the amplitude of the non- 
degenerate pole is affected by the disorder. From 
(0) and (§) we get 



J2 a NM a A^"\0) 

^E Q E k l^' k0) (k)l 2 



(18) 



Up to the terms of the second order, the sum in 
the denominator equals ^ Q \Aq%' (0)| 2 which is 
simply unity by our our choice of orthonormal 
set of Ao^(k)'s. Thus, in the non-degenerate 
case the amplitude of the pole at ko = remains 
unchanged in the first order. The amplitudes at 
other poles are: 



M (l ' k) 



where 



(19) 



should be applied. Namely, for the states with 
|k| < k c , where k c is chosen so that Sp^kc) — 

4 l) (o) « (o)Q^(o, k c ) .u_ kc A«(k c ), 

we can make the approximation Sq (k) ss s (0). 
In the first order this will not change the posi- 
tions of the poles, because k c ~ U p h and thus 
Q Q ,, ( 3(0,k c ) • Uk c ~ Uph 2 - However, all the eigen- 
vectors of the degenerate problem would contain 
Aq^(0) with finite coefficients. Thus, the unper- 
turbed weight |AfW(0)| 2 of the k = pole [see 
eq. (||)] will become distributed between all the 
poles s"' (k) that correspond to k < k c (and such 
that U_ k ^ 0). 

Consider, for example, an infinite sample in 
which each inclusion is randomly displaced from 
its original position. The sum in Eq. (^) should 
then be replaced by an integral over the first Bril- 
louin zone, so that there will be a continuous 
density of poles with nonzero weight. The cor- 
responding spectral function is written 



E4f(o)*Q Qi/3 (o,k).u_ k <j(k) 



(i,k) _ <^i_ 



F(s) 



0p 1 



1 



da', 



(21) 



(k) 



(0) 



(20) 



where |/z(s)| 2 is interpreted as a density of pole 
weight. In the purely periodic case (or for a finite 
system) , 

K( S )| 2 = ^|Af«(0)|^( S - S «(0)), (22) 



Thus, for each nonzero Fourier component Uk of 
the disorder, the poles s^ ,_k ^ in each band pick 
up nonzero amplitude M^' _k ^ ~ |Uk| ~ U p h- 

What happens in the vicinity of the unper- 
turbed poles (k = 0)? Because Q a ^(AR) — > 
as AR — > co, it follows from eq. (|12| ) that 
Q a ^(0,k — > 0) — * 0. This, however, does not 
directly relate to the small-k behavior of M^ l ' k \ 

since one also has Sq (k — > 0) — > Sq\o) at least M s ) = ^peak{s) + Hband(s), 



so that F(s) is non-analytic only at the simple 
poles s (i) (0)@. In the presence of disorder, how- 
ever, 

(23) 



for i = j. It is natural to expect that Sq (0) 
is either at the top or at the bottom of the i- 

th band, in which case Sq °(k) - 4 l) (°) ~ ak2 

for k — > 0, whereas the Q's are likely to fol- 
low Qq.flfO; k) ~ k [which would be the case if 
(|]) and ( |l3| ) were continuum Fourier transforms] . 

In this case, eq. (^0|) predicts that <4q — > oo 
as k — > 0, which means that the expression 
for the first-order non-degenerate case cannot be 
used under these conditions. Instead, for suf- 
ficiently small values of k a degenerate theory 



where n pe ak{s) represents the amplitudes of the 
discrete poles sW(k) which are degenerate with 
s^(0), while fJ.band(s) gives the amplitude den- 
sity for the rest of the band. From Eqs. ( |l9| - 
|0|), Hband{s) ~ Uph- For the typical case dis- 
cussed in the previous paragraph, the width of 
the peak part can be estimated by noting that 
s { l) (k c ) - 4°(0) - ak 2 c ~ U% h . Therefore, we ex- 
pect each 5- function entering /j,q(s) to acquire a 
half-width proportional to U? h . This typical sit- 
uation is illustrated in Fig. [l[. 
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Figure 1. Schematic of (a) the "band structure" 
for a periodic composite (two "bands" are shown) ; 
and the corresponding pole spectrum for (b) a 
purely periodic, and (c) a slightly disordered sys- 
tem (with characteristic "phononic" amplitude 
U ph ). 



dieted by the theory for the periodic case [see g 
for explicit expressions for spherical inclusions], 
but slightly broadened by a half-width propor- 
tional to T 2 ; and (b) a continuous contribution 
from a branch cut introduced by the disorder. 
The latter contribution should be most prominent 
at negative values of ei/e2 away from the original 
resonances. It would be of great interest if such a 
spectrum could be detected in a real material, e. 
g., in a suspension of metal spheres in a dielectric 
host. 
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For s sufficiently different from Sq (0), f^ pe ak{s) 
can probably be replaced by a set of <5-functions: 



F(s) 



«'6{«W(0)} 



|AfW(0)| s 



\HbandW)\' 



ds'(24) 



The second term describes the part of the spectral 
function which cannot be characterized purely by 
simple poles. Any such function can be described 
by a branch cut along the segment [0, 1) of the real 
axis, as can be seen by writing the Cauchy for- 
mula for the contour encircling the branch cut in- 
fmitesimally below the real axis. For example, the 
effective medium approximation (EMA)|| gives 
such a branch cut. While branch cut in the 
present case would certainly differ from the EMA 
one, it is quite remarkable that such a cut would 
appear even in a weakly disordered system. 

A possible application of this work could be 
a colloidal crystal at some finite temperature T. 
The total contribution to the spectral function 
from the integral in (pi|) involves \nband(s')\ 2 , 
which is determined by |U q | 2 . In a conven- 
tional crystal, the analogous sum increases lin- 
early in T at high T and approaches constant 
value as T — > 0, and similar behavior should be 
observed in a colloidal crystal. 

We thus suggest that in a colloidal crystal, at 
T such that first order corrections are adequate, 
the spectrum should consist of two parts: (a) elec- 
trostatic resonances at the same ratios ex/e^ pre- 
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